Hadronic light-by-light scattering in the muon g — 2: a Dyson-Sch winger equation approach 
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We determine the hadronic light-by-light scattering contribution to the anomalous magnetic moment of the 
muon using the framework of Dyson-Schwinger and Bethe-Salpeter equations of QCD. Our result for the pseu- 
doscalar r), rj') meson exchange diagram is commensurate with previous calculations. In our calculation of 
the quark loop contribution we improve upon previous approaches by explicitly implementing constraints due 
to gauge invariance. The impact of transverse contributions, presumably dominated by vector meson poles, are 
only estimated at this stage. As a consequence, our value a LBMuaikioo P _ ^gg _|_ 59^ x if)" 11 is significantly 
larger. Taken at face value, this then leads to a revised estimate of the total a M = 116 591 891.0(105.0) x 10" 11 . 

PACS numbers: 12.38.Lg, 13.40.Em, 13.40.Gp, 14.60.Ef 



I. INTRODUCTION 

One of the most impressive successes of the standard model 
of particle physics is the determination of the anomalous mag- 
netic moment of the electron. This quantity is determined 
both experimentally and theoretically to such a degree of pre- 
cision that the underlying physical description is vindicated. 
However, when it comes to the question of new physics, the 
anomalous magnetic moment of the muon is an even more in- 
teresting quantity, see e.g. JU-JU] for reviews. This is due to 
the large mass of the muon as compared to the electron, which 
leads to an enhanced sensitivity to physics in and beyond the 
standard model. Experimental efforts at Brookhaven and the- 
oretical efforts of the past ten years have pinned a M down to 
the 10 -11 level, leading to significant deviations between the- 
ory and experiment J3l : 

Experiment: 116 592 089.0(63.0) x 1CT 11 , (1) 
Theory: 116 591 790.0(64.6) x 10~ u . (2) 

Whilst the theoretical and experimental values are determined 
to comparable errors, the central values give rise to a discrep- 
ancy at the 3.3cr confidence level. This difference has been 
present for a number of years and can be interpreted as a signal 
for the existence of physics beyond the standard model. How- 
ever, to clearly distinguish between New Physics and possible 
shortcomings in the SM calculations the uncertainties present 
in both experimental and theoretical values of a M need to be 
further reduced. 

The greatest uncertainties in the theoretical determination 
of a M are encountered in the hadronic contributions, i.e. those 
terms which involve QCD beyond perturbation theory. The 
most prominent of these is given by the vacuum polarisation 
tensor dressing of the QED vertex, see Fig. |l(a)| Fortunately 
it can be related to experimental data of e + e~ -annihilation 
and r-decay via dispersion relations and the optical theorem, 
thus resulting in a precise determination with systematically 
improvable errors Q|]. Considered individually, its (leading 
and subleading order) contribution to the anomalous magnetic 
moment of the muon is JH] 



FIG. 1: The two classifications of corrections to the photon-muon 
vertex function: (a) hadronic vacuum polarisation contribution to a^. 
The vertex is dressed by the vacuum polarisation tensor I1 M „ ; (b) the 
hadronic light-by-light scattering contribution to a M . 



Although currently these uncertainties dominate the error of 
the theoretical result in Eq. (|2]i it is foreseeable that future ex- 
periments reduce this error below that of another, more prob- 
lematic source. This is the hadronic light-by-light (LBL) scat- 
tering diagram, shown in Fig. |l(b)| This contribution cannot 
be directly related to experiment and must hence be calculated 
entirely through theory. The central object in such a calcu- 
lation is the photon four-point function. It receives important 
contributions from the small momentum region below 2 GeV, 
where perturbative QCD breaks down and non-perturbative 
methods are imperative. Recent determinations of aJ^ BL are 
provided in Table U Although the magnitude of the LBL con- 
tribution is much smaller than the one from vacuum polar- 
ization it is significant because its error is of a comparable 
size. Taken together, with the errors added in quadrature, the 
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TABLE I: Recent calculations of the hadronic light-by-light scatter- 
ing contribution to the anomalous magnetic moment of the muon. 
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hadronic contributions constitute the largest uncertainty in the 
standard model determination of the anomalous magnetic mo- 
ment of the muon. 

The theoretical approaches to determine the LBL contribu- 
tion are centered around two main ideas. One is chiral sym- 
metry, its breaking pattern and the associated low energy ef- 
fective descriptions of QCD flUl : the other is the large- N c 
expansion of the four-photon function and the associated or- 
dering of diagrams. These ideas have been put together in 
|[l2ll and led to various refined calculations of LBL within the 
frameworks of large- N c and vector meson dominance H, [lj- 
Il5ll . the extended Nambu-Jona-Lasinio model (ENJL) i fflflql . 
the (very similar) hidden local symmetry modeljl8T- [l0ll . or a 
non-local chiral quark model lfl7[[l8ll . see also Jajfor a sum- 
mary. Although in terms of diagrams individual contribution 
are of varying size in these approaches, their sum leads to con- 
sistent results as can be inferred from Table |U In all these cal- 
culations the (pseudoscalar) meson exchange contributes the 
most and the meson loop has been found to be small. A pos- 
sible explanation of the latter is given in Hill . As a result, 
we quote the recent value for LBL a^ BL = 105(26) x 10~ n 
proposed in Ref. which also agrees with the one in 

One of the most important goals for these and future cal- 
culations is the reduction of the model dependence and sub- 
sequently of the systematic error involved in these calcula- 
tions. Since LBL is non-perturbative in nature, all estimates 
in Table U are plagued by systematic model dependencies. It 
is therefore desirable to also explore other calculational tools 
which have the potential to go beyond these limitations. Cer- 
tainly, lattice gauge theory is one such method. However, 
due to the multi-scale nature of the problem no reliable es- 
timates for LBL have been extracted on the lattice so far. This 
multi-scale nature also makes EFT methods less desirable as it 
proves more difficult to impose suitable matching conditions. 

Another non-perturbative method, well suited to accommo- 
date for largely different scales is the framework of Dyson- 
Schwinger and Bethe-Salpeter equations lfl9l - l22Tl . In the past 
years this approach has been used to study fundamental prop- 
erties of QCD such as confinement and dynamical chiral sym- 
metry breaking. On the other hand the approach served as a 
tool for hadron physics. In this work we expand upon this 
and apply the formalism to a calculation of the LBL contri- 
bution to the muon anomalous magnetic moment. To this end 
we separate different contributions to the light-by-light four- 
point function according to their topology of gluon exchange 
and their status with respect to the large-A^ expansion. Dia- 
grammatically this translates to considering resummations of 
planar diagrams involving gluon exchange. In this scheme 
we then determine the dressed quark-loop diagram and an 
approximation in terms of pseudoscalar meson (ir°, r], 7/) ex- 
change contributions. In principle, the off-shell meson ampli- 
tudes involved in these diagrams could be calculated from in- 
homogeneous Bethe-Salpeter equations. Here, due to numer- 
ical complexity, in this work we resort to a commonly used 
ansatz that extrapolates on-shell wave functions. Our results 
are then compared with the ones of previous approaches. First 
results of our analysis have been published in Ref. l23ll . Here 
we discuss our method in much more detail and present new 



and more elaborate results for the quark-loop diagram. 

The paper is organized thus: in section [II] we recall the 
definition of the light-by-light scattering amplitude and fo- 
cus upon the pseudoscalar pole contributions; in section [III] 
we introduce our Dyson-Schwinger approach and discuss the 
necessary truncation schemes; in section [IV] we present and 
discuss our results. We conclude in section [Vl 



II. THE LBL SCATTERING AMPLITUDE 

In the hadronic light-by-light scattering contribution, 
Fig. |l(b)[ the muon is coupled to an external photon source 
via the hadronic photon four-point function II M „ a) 3, defined 
through 



a iqi-x+iq 2 -y+iqa-z 



{j^)jv{x)j a {y)jp{z)) 



(4) 



where j xyz = J d 4 x J d 4 y Jd 4 z represents integration over 
four-dimensional space, 91,2,3 are the photon momenta that 
are connected to the muon line, and is the electromagnetic 
quark current 
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A detailed discussion of this object can be found in the liter- 
ature, see e.g. |[l4i[l6ll . Instead of working directly with the 
light-by-light scattering diagram given in Fig. |l(b)[ it is more 
convenient to follow the strategy employed in Ref. f24l l25ll . 
Here gauge symmetry is exploited to construct quantities 
that are finite. Through use of the Ward-Takahashi-identity 
k^Tl^ap = it follows via differentiation that 



n 



k — n 



(6) 



which serves as definition of the five-point-function n^^j,^. 
Here k = qi + q 2 + is the momentum of the external pho- 
ton. The virtue of the derivative is that it lowers the dimen- 
sionality of the integral thus ensuring that integrals employing 



n 



are manifestly convergent. We define the quantity 
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which is now related to the dressed muon vertex of 



Fig. [1(b)] via 



(8) 
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FIG. 2: The hadronic light-by-light (LBL) scattering contribution to a M and its expansion, using EFT approaches, as a quark loop part (left), 
leading pseudoscalar meson exchange part (middle) and a leading meson loop part (right). Note that the quarks here may be interpreted 
differently to those in Fig. [3] 
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FIG. 3: The hadronic light-by-light (LBL) scattering contribution to a M and its expansion, using functional methods, as a quark loop part (left), 
a ladder exchange part (middle) and a ladder ring part (right). All propagators and vertices are fully-dressed, with the ellipsis marks indicating 
that an infinite number of gluons are resummed. 



Here D^(q) are perturbative photon propagators (we use 
Feynman gauge) with momenta <jj. The perturbative muon 
propagators are given by S(p). 

The anomalous magnetic moment can now be obtained by 
applying the appropriate projection operator to Eq. (|7]i 
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that we write here in Euclidean convention for later conve- 
nience. Using Eqs. (((JO we are able to evaluate the light- 
by-light scattering contributions for an arbitrary photon four- 
point function. What remains now is the specification of this 
four-point function within our approach. 



A. Expansion using EFT approaches 

As already mentioned in the introduction, chiral and large- 
N c arguments have been established to expand the full LBL 
scattering amplitude into the diagrammatic parts shown in 
Fig- E] lit]. These diagrams belong to different orders with 
respect to chiral and large- N c counting. Whereas the meson 
exchange diagrams and the quark-loop diagram are leading in 
large-A^, it is the meson-loop diagram that is leading in the 
chiral counting. Thus a priori, one does not know which ex- 
pansion is to be preferred. Therefore it is certainly interesting 
that all explicit calculations of these contributions seem to fa- 
vor the A^-counting scheme; meson-loop contributions have 
been found to be suppressed. Arguments as to why this is the 
case have been presented in Ref. fl4tl . 

Strictly speaking, however, one does not actually perform 
a large-A^ expansion as this would necessitate the inclusion 



of an infinite number of resonances. Instead, only the low- 
est lying meson exchange contributions in the pseudoscalar, 
scalar and axialvector channel have been subsummed. Here, 
the pseudoscalar 7r°-exchange has been identified as the lead- 
ing contribution, followed by rj and ^'-exchange. 

Concerning the pseudoscalar (PS) exchange contribution a 
few remarks are in place. The photons in the exchange di- 
agrams are coupled to the PS mesons via the PS-77 form- 
factor, fps77- It is evident that there are two limiting features 
of the pseudoscalar-pole approximation. The first is the actual 
provision of the form-factors themselves, which are in general 
subject to systematic errors depending on how they are mod- 
eled or calculated. The second is the procedure under which 
the form-factor is taken off-shell. Previous approaches mainly 
used vector meson dominance ideas to determine this form 
factor and there has been an extensive debate as to whether 
and how short distance constraints have to be implemented 
III 111 Q31 [Still]. Rather than employ the principles of vector- 
meson dominance and construct an ansatz for the on-shell/off- 
shell form-factor, we wish to calculate it from first principles. 
This is possible within the framework of Dyson-Schwinger 
and Bethe-Salpeter equations using a well explored and suc- 
cessful truncation scheme 1201. 



As for the quark-loop diagram, different interpretations 
have been given in the literature. Whereas in lfl2ll it has been 
argued that the quark loop is a separate contribution that has 
to be added to the other two, in many other approaches it has 
been treated as a complementary one, which is only added 
in the large spacelike momentum region, say above a typi- 
cal cutoff for an effective model. In our functional approach, 
described below, it is clearly the first point of view that is cor- 
rect. Moreover, as we will see, the quark loop is subject to 
large dressing effects not only for the quark propagators in the 
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loop but also for the quark-photon vertices. This will be the 
main result of our work. 



B. Expansion using functional methods 

From a functional integral approach to QCD, featuring 
quarks and gluons as the fundamental degrees of freedom, the 
analogous picture to what is normally considered in the lit- 
erature is shown in Fig. [3] where we give an expansion in 
terms of non-perturbatively dressed one-particle irreducible 
Green's functions. The basic idea of this expansion is not 
a separation of long distance and short distance scales, but 
rather a separation of different classes of diagrams based on 
their topology. Clearly, the expansion is such that no dou- 
ble counting of diagrams is involved. By considering a re- 
stricted subset of contributions in which only diagrams with 
a planar topology are resummed, we effectively adhere to the 
A^-counting scheme as favored in the EFT approaches men- 
tioned above. Though there are similarities between the two 
pictures, since we work with a truncated formulation of exact- 
QCD rather than an effective field theory there are some dif- 
ferences that we will comment on here to avoid confusion. 
First of all, our quarks are to be interpreted in the same way 
as those extracted via Lattice QCD; they are characterised by 
momentum dependent dressing functions that interpolate be- 
tween the current and constituent quark limits, cf. the discus- 
sion below Fig.QT] Secondly, the quark-photon coupling is a 
non-perturbative form-factor and not merely a tree-level bare 
vertex; it can be calculated self-consistently for a given trun- 
cation scheme. Finally, the planar resummation of gluons is 
related to the T-matrix of quark-antiquark scattering and con- 
tains meson poles that can be associated with pseudoscalars, 
vectors, scalars etc. This will be exploited below, where we 
return to the conventional meson exchange picture to approx- 
imate these contributions. 

We wish to emphasize that the expansion displayed in Fig. [3] 
has been used successfully in a different context already in 
Ref . ll26ll . There it — ir scattering has been considered using 
similar quark-box and ladder exchange part s as displayed in 
Fig. [3] In this setup, the authors of Ref. 1 26 1 could reproduce 
the isospin and 2 scattering lengths in exact agreement with 
Weinberg's low energy results. Moreover, in Ref. ll27ll it has 
been checked, that the corresponding resonant expansion sim- 
ilar to the one displayed in Fig.|2]is a good approximation to 
the ladder exchange part of Fig. [3] Note that in both these cal- 
culations the quark-box diagram had to be added to the ladder 
exchange or the resonant 'meson-exchange' part respectively. 
We believe that these results add further support to our ap- 
proach. 



1. Quark-loop contribution 

Within our proposed truncation, the quark-loop is com- 
posed of dressed quark propagators and dressed quark-photon 
vertices. On expanding these one-particle irreducible Green's 
functions, within the rainbow-ladder approximation, we find 




FIG. 4: Expansion of quark-loop contribution to the photon four- 
point function in terms of planar quark and gluon diagrams (all prop- 
agators are fully dressed). 




FIG. 5: The hadronic light-by-light scattering contributions to a,, 
from the quark loop. There are an additional three diagrams (not 
shown) in which the quark spin-line is reversed. Principally, these 
diagrams involve dressed quark propagators and quark-photon ver- 
tices. 



planar-like diagrams such as the ones shown in Fig. [4] (all 
propagators are fully dressed), where in fact infinite ladders 
of gluons are taken into account. Should we consider cor- 
rections beyond rainbow-ladder, such as those considered in 
Ref. l28l l29ll . one would also include diagrams in which the 
gluons have self-interactions as well as crossed-ladder com- 
ponents. Taking into account such corrections is, however, 
beyond the scope of the present work. 

Considering this contribution to the muon-photon vertex, 
we obtain the diagrams as shown in Fig. |5j where we have 
shown permutation of the external photon legs but have omit- 
ted the topologies that merely involve reversal of the quark- 
spin line (these give identical contributions and hence con- 
stitute a factor of two). As is well-known, these diagrams 
are individually logarithmically divergent with only their sum 
finite and convergent; thus one employs the aforementioned 
trick, Eq. @, of taking the derivative of the photon four-point 
function I24II25I1 

Since this is now to be applied to loop integrals over 
non-perturbative quantities, namely the quark-propagator, it 
is no longer possible to reduce the integration to be five- 
dimensional as in the case of perturbative studies. More 
generally, on considering the planar nature of the diagrams, 
one must deal with 8-dimensional integrals which necessi- 
tate Monte-Carlo methods 13011 . However, for reasons of cal- 
culational simplicity we actually integrate in nine. We did 
check, however, that we were able to reproduce the well- 
known perturbative results for the electron loop contribution 
to the anomalous magnetic of the electron and the muon ll3~ll - 
l33ll . Additionally, due to the somewhat involved Dirac alge- 
bra ll34l l35ll we will content ourselves with taking the quark- 
photon vertices inside the quark-loop contribution to be: (a) 
bare; (b) 1BC; (c) full BC. The precise meaning of these ab- 





FIG. 6: Ladder-exchange contribution (upper equation) and ring- 
ladder contribution (lower equation) to the photon four-point ampli- 
tude. 




FIG. 7: The T-matrix in Rainbow-Ladder approximation, (a) shows 
the series expansion in terms of dressed quarks and gluon, whilst 
(b) represents Dyson's equation, (c) Shows the pole-ansatz for the 
T-matrix on-mass shell. 



breviations and the relation to the full quark-photon vertices 
will become clear in sections Ull B I and \l V B I The extension to 
employ the numerically calculated non-perturbative form of 
the vertex will be explored in a later publication. The results 
of our calculation are presented in section ITVl 



2. Ladder-exchange and ladder-ring contribution 

Two contributions that are leading and sub-leading in large- 
N c respectively are the so-called ladder-exchange and ladder- 
ring diagrams of Fig. [3] These infinite ladder resummations 
are in fact related to the T-matrix of bound-state theory in 
a certain approximation scheme (that produces planar dia- 
grams). Thus, another way to portray these contributions is 
given in Fig. [6] The T-matrix in rainbow-ladder approxima- 
tion is given in Fig. [7] At this point, we make it clear that there 
is no conflict nor double counting between the quark-loop and 
ladder-exchange diagrams, as they clearly consider and resum 



FIG. 8: Pole representation of the ladder-exchange contribution to 
the photon four-point function. 




FIG. 9: The pion-pole part of the LBL contribution to a M . The three 
possible permutations of the photon legs are not shown. 

different topologies of diagrams. 

As it stands, the full T-matrix is a very complicated object 
to solve in its entirety though its structure admits several ap- 
proximations and simplifications l36ll . The one which we em- 
ploy here is similar to the viewpoint taken by Effective Field 
Theory approaches; that is, we consider pole contributions to 
be dominant. Now, since it is well-known that such an infinite 
gluon-ladderresummation dynamically generates bound-state 
poles, one can expand the T-matrix in terms of meson pole 
contributions as shown in Fig. |7jc). On mass shell we then 
have a unique definition of the Bethe-Salpeter amplitude, de- 
scribed below in Fig.[l2l that gives the form-factor describing 
coupling of a meson to two quarks. From this point the (on- 
shell) pseudoscalar-photon-photon form-factor can be defined 
and calculated, giving rise to the 'leading' pseudoscalar me- 
son exchange part, as shown in Fig.[8]and Fig. [9] 

In a similar fashion, the ring-ladder diagram contains con- 
tributions akin to the pion loop on meson mass-shell. How- 
ever, since these are generally considered to be sub-leading 
we will not consider them further here and instead concen- 
trate on the quark-loop and ladder-exchange diagrams. After 
we present our approach and formalism in the next section, 
results will be discussed in section HVl 



III. FRAMEWORK 

The dressed quark propagator is one of the most impor- 
tant quantities in the covariant description of mesons. It en- 
codes non-perturbative properties of QCD such as dynamical 
mass generation and the realization of a non-zero condensate. 
Its equation of motion, the quark Dyson-Schwinger equation 
(DSE) displayed in Fig. [TO] also contains the dressed gluon 
propagator and a dressed quark-gluon vertex. Whereas the 
dressed gluon propagator in Landau gauge is a well-known 
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FIG. 10: Dyson-Schwinger equation for the quark propagator. Spec- 
ification of the fully-dressed gluon propagator (wiggley line) and 
quark-gluon vertex (coloured blob) defines the truncation scheme. 



quantity by now, see IB7H40I1 and references therein 1 , the 
study of the details of the dressed quark-gluon vertex is still on 
an exploratory level, although some progress has been made 
in the past years ll29i l42l 14311 . Pending deeper insights into 
the nonperturbative structure of the quark-gluon interaction it 
is therefore reasonable to work with approximations that take 
into account important features of the full theory. This strat- 
egy, of course, introduces model dependencies into our calcu- 
lation that have to be carefully addressed later on. 

>From a phenomenological perspective a successful ap- 
proximation in this respect is the rainbow-ladder truncation of 
the quark-DSE. The philosophy here is to combine the dress- 
ing of the gluon propagator with the vector part of the quark- 
gluon vertex into a single function depending on the gluon 
momentum only. While this is certainly a severe approxima- 
tion in principle, in practice it turned out to be very success- 
ful as concerns the calculation of masses and electromagnetic 
properties of mesonic observables 12011 . While the parame- 
ters of the model are tuned such that it reproduces the experi- 
mental values for the masses and decay constants of the pion, 
it also reproduces the pion charge radius and 7177 transition 
form factors on the percent level. In the vector channel the 
agreement with experimental masses and decay constants is 
on the five and ten percent level. Thus, while one has to keep 
in mind possible systematic caveats, we nevertheless believe 
that such a model is an excellent starting point for a systematic 
evaluation of hadronic LBL. 

In Euclidean momentum space, the renormalized dressed 
gluon and quark propagators in the Landau gauge are given 
by 



D^uip) = ( <V 2" 



Sf(j>) 



j p* 



(12) 



iff + M(p 2 ) iffA(p 2 ; p 2 ) + B(p 2 \ /I 2 ) ' 



(13) 



where Z(p 2 ;/i 2 ) is the gluon dressing function, Zf(p 2 ;/j, 2 ) 
is the quark wave-function and M(p 2 ) is the renormalisation 
point independent quark mass function. The dependence of 
such functions on the renormalisation point pi 2 will be im- 
plicitly assumed from here on. The quark dressing functions 



1 There is an intense debate on the behaviour of the gluon propagator in the 
deep infrared, i.e. for momenta p < 50 MeV. It seems, however, that 
this momentum region is irrelevant when it comes to the calculation of 
observables 128,29,411. 
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FIG. 11: Mass function and wave-function dressing functions corre- 
sponding to quark propagators solved with the Maris-Tandy interac- 
tion |3. 



A{p 2 ) and B(p 2 ) can be recombined into the quark mass and 
wave-function by M{p 2 ) = B(p 2 )/A(p 2 ) and Zf(p 2 ) = 
1/A(p 2 ). 

These propagators may be obtained by solving their re- 
spective Dyson-Schwinger equations. The DSE for the quark 
propagator, shown diagrammatically in Fig.[l0] is written 



S- 1 (p) = Z 2 So 1 (p) + ^( P ) , 

—^T u (q,p)D^(k)^S F (q) , 



(14) 



where is the quark self-energy, k = p—q and the Casimir 
Cf = 4/3 stems from the colour trace. We introduced the 
reduced quark-gluon vertex T v (q,p) defined by T®(q,p) — 
ig^-T v {q,p). The bare inverse quark propagator is Sq 1 (p) = 
iff + 771. The renormalisation factors are Z\p = Z2/ Z3 for 
the quark-gluon vertex, Z2 for the quark propagator and Z3 
for the ghost dressing function. 

The scalar dressing functions of the quark DSE are solved 
for by appropriate projections of Eq. ( TBI ). This is a coupled 
non-linear integral equation that is solvable provided we know 
the gluon dressing function and the structure of the quark- 
gluon vertex. In the rainbow approximation both are spec- 
ified by Ansatze, with in particular the choice T v (q,p) := 
T YM (k 2 )j u , with scalar function r YM representing the non- 
perturbative dressing of the quark-gluon vertex and k = q — p 
the gluon momenta. Here, the gluon dressing function Z(k 2 ) 
from Eq. (Q2I and the Yang-Mills part T YM (fc 2 ) of the quark- 
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gluon vertex are combined to form a phenomenological effec- 
tive interaction. For the Maris-Tandy (MT) model i44ll this 
function is given by 

z(fc 2 )r YM (fc 2 ) = ^^m 4 cx P (-fc 2 /^ 2 ) 

27T7 m 

log (r+(l + fcVA 2 QCD ) 2 ) 

x [l-cxp(-fc 2 /[4m 2 ])]) ,(15) 

with 

mt = 0.5 GcV , r = c 2 - 1 

7m = 12/(33 -27V/) , Aqcd = 0.234 GcV. 

This interaction corresponds to a Gaussian distribution in the 
infrared that provides for sufficient interaction strength to gen- 
erate DCSB, together with the one-loop behavior of the run- 
ning coupling at large, perturbative, momenta. The latter is 
mandatory to provide for the correct short distance behavior 
of the quark propagator. The remaining parameters u and 
D essentially constitute a single one-parameter family of so- 
lutions for which pion observables remain comparable, via 
ujD = (0.72GeV) 3 . 



For the convenience of the reader, in Fig.QT]we again show 
the two dressing functions Zf(p 2 ) and M(p 2 ) that character- 
ize the non-perturbative quark propagator, obtained by solv- 
ing Eq. ( f]~4T > using the Maris-Tandy interaction l44ll . Clearly, 
in the mass function there are three distinguished momen- 
tum regions. In the infrared, the quark propagator is essen- 
tially constant displaying the behavior of a constituent quark. 
Then for 1 GeV 2 < p 2 < 10 GeV 2 there is a region of 
rapid change, where the quark mass function follows the well 
known 1/p 2 behavior expected from the operator product ex- 
pansion. For even larger momenta and non-vanishing current 
quark mass the quark mass function behaves logarithmically 
as expected for a current quark. The fully dressed quark prop- 
agator thus naturally interpolates between the constituent and 
current quark picture. We consider this feature of the Dyson- 
Schwinger approach to QCD as an advantage compared with 
effective models such as the ENJL model. 



A. Bethe-Salpeter equation 

The chiral symmetry preserving truncation for the Bethe- 
Salpeter equation, consistent with the rainbow-approximation 
above, is given by the ladder approximation 



K tu . rs ( P , k; P) [S F (k+)T"(k; P)S F (k-)] 



with the kernel K tu . rs given by 



g 2 Z{k 2 )T™(k 2 )Z lF ( k^ 
K tu . }Sr (q,p;P) = -J IV p- 











) 




ts 


T 7v . 



(16) 



(17) 



see Fig. Q~2] for a graphical representation. Here T qq (p; P) is 
the Bethe-Salpeter vertex function corresponding to a pseu- 
doscalar quark anti-quark bound-state, specified below. The 
momenta k + = k + P/2 and fc_ = k — P/2 are such that the 
total momentum P of the meson is given by P = k + — fc_ 
and the relative momentum k = (k + + k—)/2. The Latin 
indices (t, u, r, s) of the kernels refer to colour, flavour and 
Dirac structure. 

The form of the kernel Eq. ( flTT ) is uniquely determined from 
the axial-vector Ward-Takahashi identity and ensures, that the 
pion is a Goldstone-boson in the chiral limit without any fine- 
tuning of parameters. It also ensures that important constraints 
from chiral symmetry such as the Gell-Mann-Oakes-Renner 
relation are satisfied. 



doscalar meson is completely specified by the following form 



T™(p; P) = 75 |fi w (p; P) - W(p; P) 

~irf{p- P) Ff{p; P) - \P^] Ff(p; P) 



(18) 



This amplitude is obtained through solution of Eq. ( TToT ) on- 
mass shell: P 2 = —m 2 q in Euclidean space. While (TT~8b rep- 




In general, the covariant structure of the Bethe-Salpeter ver- 
tex function, T q ^{p\ P), determines the quantum numbers of 
the bound-state under consideration. In particular, a pseu- 



FIG. 12: The homogeneous Bethe-Salpeter equation for the meson 
amplitude. 
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FIG. 13: Normalised Bethe-Salpeter amplitudes Fi and F2 of the 
pion solved with the Maris-Tandy interaction (leading Chebyshev 
component) 14411. 



resents a quark-anti-quark bound-state, physical mesons are 
defined as matrices in flavour space built out of the qq am- 
plitudes. This then leads to the same decomposition as in 
Eq. ( fT8l but with flavour matrix valued quantities V and In 
the following, however, we will keep the flavour index implicit 
and use Eq. (fT8b for qq amplitudes and mesons alike. Anyway, 
in the isospin-limit considered herein the pion amplitude dif- 
fers from the uu/dd amplitudes only by flavour-matrix struc- 
ture. The pole masses and the scalar amplitudes that contain 
the dynamical information are identical (up to normalization). 
This is different for the 77 and 17'. 

In the chiral limit the leading behavior of the pion amplitude 
T 77 is given by 



(19) 



where B(p 2 ) is the scalar dressing function of the quark, f„ 
is the chiral limit value of the leptonic decay constant, and 
A3 is a Gell-Mann matrix that represents the flavour structure. 
The hat in Eq. (fT9] l indicates that the object is matrix valued in 
flavour space. The expressions for the calculation of f„ and 
the normalization condition of the Bethe-Salpeter amplitude 
together with details on the numerical procedure for dealing 
with the BSE are given in ll45ll . For the convenience of the 
reader we display the resulting Bethe-Salpeter amplitudes for 
p 2 = —in 2 and p.P = in Fig.[T3] Qualitatively, the am- 
plitudes have a similar form as the quark mass function in the 
chiral limit. For large momenta and up to logarithmic correc- 
tions they fall off like l/p 2 , which is a necessary condition 
to correctly describe the anomalous decay of the pion and to 
reproduce the asymptotics of the pion form factor [46]. 



B. Quark-photon Vertex 

An important quantity for the determination of the LBL 
contribution to the muon g — 2 is the dressed quark-photon 
vertex. This quantity is genuinely non-perturbative in nature 
and necessary for the calculation of the PS77 form-factor. It 



describes the coupling of a fully dressed quark to a photon 
and is dominated by QCD corrections. As a function of one 
Lorentz and two spinor indices, it can be decomposed into 
twelve Dirac structures 



(20) 



where V*(P, k) represents the basis components, and \i(P, k) 
the non-perturbative dressing functions. A common basis is 
that of Ball and Chiu j47ll . in which VHP, k) is split into terms 
that are transverse and non-transverse with respect to the pho- 
ton momentum. The Ward-Takahashi identity and regularity 
assumptions constrain the form of the non-transverse part in 
terms of quark propagator functions 



C 3 := 



7m" 



A{k 2 ) + A{k 2 _) 



1 A(k 2 
2 



A(k 2 _) 



k\ 



~ l (fc+ + fc_) i 



B{k\) - B(k 2 



k 2 _ - k\ 



(21) 



leaving only the strictly transverse pieces undetermined. 
Eq. ( t2Tb can therefore be seen as an approximation of the full 
vertex and has been used in situations, where the full vertex 
cannot be determined numerically. 

A more sophisticated approach, however, is to solve the 
inhomogeneous Bethe-Salpeter equation of the quark-photon 
vertex, shown diagrammatically in Fig. [14] It is given by 

Tt u {k,P) = Z ll L + J K tu . irs [S F (k + )T^q,P)S F (k.)] sr 

(22) 

where Z\ is the renormalisation factor associated with the 
quark-photon vertex. By using the same interaction kernel as 
in the BSE for mesons, Eq. (1171) . not only do we achieve self- 
consistency within the truncation scheme, but also by virtue of 
its symmetry preserving nature we satisfy the Ward-Takahashi 
identity. Consequently, the non-transverse part of the vertex, 
given in Eq. (121 l l. is nicely reproduced numerically, with trans- 
verse terms additionally generated ll48l l52U54ll . 

We note here that such a determination of the quark-photon 
vertex automatically contains poles in the time-like region 
corresponding to vector meson exchange. Thus, presuppos- 
ing that vector-meson dominance is an important feature in 




+ 




FIG. 14: Inhomogeneous BS equation for the quark-photon vertex, 
in rainbow-ladder approximation. 
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the structure of the pion electromagnetic form-factor, it is al- 
ready included here as a result of the approach we employ. 
This has been discussed in detail also in Refs. HHHD]. 

The numerical details involved in the calculation of the 
quark-photon vertex have been described in several works, see 
e.g. the appendix of Ref. HH . Below we will use the fully dy- 
namical, selfconsistent solution of Eq. d22l ) for our calculation 
of the PS — > 77 form-factor and the resulting meson exchange 
contribution to LBL. Unfortunately, because of the numerical 
complexity we have to restrict ourselves to the exact longitu- 
dinal part given by Eq. ( f2TT > in the quark-loop diagram. It will 
be a subject of future work to overcome this limitation. 



C. The PS — ► 77 form-factor 

The coupling of the exchanged pseudoscalar mesons to 
photons is the quantity that is central to the resonant expan- 
sion of Fig. [2] In impulse approximation, consistent with the 
rainbow-ladder truncation scheme introduced in sections [Till 
the Bethe-Salpeter amplitude is connected via a quark-triangle 
to the fully-dressed quark-photon vertex, as shown in Fig. [15] 
For a pseudoscalar PS we have 

A^'->'{k 1 ,k 2 ) =2e 2 N c J \i[iQ e T4k 2 , Pl2 )S F ip 2 ) 

k 

xf PS (p 23 ,P)S F (p 3 )iQ e T, l (k u p 31 )S F ( P3 )] , (23) 

where k% and k 2 are the outgoing photon momenta, p\ = q, 
p 2 — q—k 2 and 753 = q + ki are the quark momenta and pij = 
(pi + Pj) /2. The factor of two stems from exchange of the 
two photon vertices and Q = diag [2/3, —1/3,-1/3] gives 
the quark's charge. The PS vertex T PS is explicitly matrix 
valued in flavour space. It is defined as 



1 

i 

7! 



diag 



diag 



pirn pdtZ q 

j^uu pcid op 
pirn pdeZ p 



(24) 



for the pseudoscalar mesons. The T qq are solutions of 
Eq. ( TToT l. In addition we work in the isospin-limit {T uu = 
T dd ). Since the quantities in Eq. (l24l) are defined in the singlet- 
octet basis we have to rotate in order to obtain the 77-77' ampli- 
tudes 



f =cos6>r" -sin 9 T" 



L0F> 



(25) 



where we haven taken 9 = —15.4° lf55ll . The pseudoscalar 
electromagnetic form-factor can be described by a single 
scalar function, F" 7 7 . For the pion this function can be 
given a natural normalization via the Abelian anomaly j5<l] 



(26) 




FIG. 15: The 7r°77 form-factor in impulse approximation. All inter- 
nal quantities are fully-dressed. 




FIG. 16: The ir — > 77 form-factor plotted as a function of the two 
photon momenta k\ and We compare our numerical results in the 
Dyson-Schwinger/Bethe-Salpeter approach with an ansatz inspired 
by vector meson dominance discussed in Ref. [5J]. 



where a cm is the fine structure constant and f n the pion de- 
cay constant. The definition of the prefactors is such that 
^77(0,0) = 1. The 77- and ?7'-mesons have the same ten- 
sor structure. 

Note that the form factors determined here do not accu- 
rately reflect all effects due to the topological mass of the 
770, simply because the J7 y i(l)-anomaly is not represented cor- 
rectly in the Maris-Tandy model 2 . In the form factors this 
may be a minor problem. The effect is larger, however, in 
the meson propagators attached to the form factors. We there- 
fore prefer to use the experimental masses in these propaga- 
tors thereby taking care of the majority of the [/^(^-anomaly 
effects. 

The 7T electromagnetic form-factor has been explored in 
detail in Ref. ll53ll . wherein it has been confirmed that the cor- 
rect normalization is satisfied. In addition it has been shown 
analytically (and numerically) that the correct asymptotic be- 
haviour, modulo potential logarithms is obtained f53il57ll . 



lim F 

Q 2 ->oo 



lim 7 * 7 *(Q 2 ,<2 2 ) oc — 

Q 2 -s-oo 



(27) 



In Fig [16] we plot the form factor as a function of the two pho- 
ton momenta kf and k\ and compare with the VMD inspired 



2 Perspectives to improve this issue in the framework of Dyson-Schwinger 
equations have been reported in Refs. l4^ - [5lll . In Ref. t5lll a topological 
mass of the r]o has been obtained which goes well with lattice results of the 
topological susceptibility via the Witten-Veneziano relation. 
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model used in Ref. 11311 . One clearly sees that both form fac- 
tors agree nicely on a qualitative and even quantitative level. 
Whereas the low-momentum behaviour is governed by the 
anomaly, at large momenta both form factors fall off accord- 
ing to Eqs. ( f27l >. There are small quantitative differences in 
the mid-momentum regime, which will lead to a small differ- 
ence in the meson exchange contributions to LBL, discussed 
below. In general, however, the results of our calculation may 
be viewed as a confirmation of the previously used model ap- 
proaches almost from first principles. 



Here k± = k ± P/2. Note that in the chiral on-shell limit 
(P 2 = 0) the above equation reduces to Eq. 1 1191 . We gen- 
eralize the pseudoscalar amplitude Eq. <fl~8b by using Eq. d29l ) 
as a guideline for all four structures also away from the chiral 
limit. The final off-shell meson amplitude reads: 



1 =75 



pPS 



(p;P) + f(P 2 ){ 



itf(p-P)FP{p-P)-\PJ]Fl s 



(30) 



D. Off-shell prescription 

It is evident from the kinematics of the diagram shown in 
Fig.|9]that the form-factors, thus far defined as on-shell quan- 
tities, must be evaluated for momenta of the exchanged pseu- 
doscalar meson that would be far from the pole mass. In 
the approach considered here, the pseudoscalar amplitude is 
obtained from its homogeneous Bethe-Salpeter equation and 
hence is by definition an on-shell quantity. Thus, to proceed 
we must introduce a prescription for the continuation of this 
quantity to the off-shell momentum region. 

Since the off-shell behaviour should be dominated by the 
pseudoscalar pole contribution, the introduction of any pre- 
scription that provides for a suppression at off-shell momen- 
tum should be suitable starting point. Here, we will employ 
a prescription that is inspired from the axial-vector Ward- 
Takahashi identity in the chiral limit 

2P^ 3 (k,P)=iS- 1 (k + ) l5 + i l5 S- 1 (k_) . (28) 

Here the axial vector vertex is defined as the correlation func- 
tion r^' 3 = (Jn' 3 QQ) that includes the axial vector current in 

the pion channel j^ 3 = qj^"/ 5 ^-q. It is clear that Eq. ( |28l 
relates T^' 3 to the quark propagator. Taking explicit param- 
eterisations for vertex and propagator (see 114511 ) that include 
the pion pole in the axial vector vertex, the following form of 
the dominant amplitude for the tt° can be deduced: 

ff(t ,P,. A! *l±ffi, (29) 



where the hat over the functions Ff s indicates that the flavour 
structure of the corresponding meson is included in the same 
manner as in Eqs. d24l |25| ). The scalar off-shell amplitudes 
F 99 (k,P) are defined in terms of the on-shell amplitudes 3 
Ff (k, k ■ P) through 

^,P) = ^^' f 'f( L ' L ' f ) (31) 

for which i = 1, . . . , 4. The on-shell amplitudes are obtained 
via Eq. (TToT l. In order to account for the mass dimensions of 
the form-factors i^ff 4 we attach to each the function 

This prevents an unnatural enhancement at high meson virtu- 
ality whilst at the same time leaving the on-shell behaviour 
unchanged. The off-shell form factor 

Al^'(P,h,k 2 ) = eM ^fc«^ PS . 7 . r (P 2 ,fc2 ; fc 2 2 ) , 

(33) 

is then obtained via the generalisation of Eq. d23l ) by taking 
the Bethe-Salpeter amplitude T 9 ^ to be defined via Eq. ( l30l . 

The contribution to the derivative of the four-point-function 
can now be written as II 1311 : 



fr / s •7 r PS* 7 *7* (9l2:9l:92)-^PS*7*7*(9l2:<7l2> ) a p T 

ii (p)^Ao-(9lj 92, — qi2) — 2~T 2 ^nvafill <?2 e Ao-pr9l2 

q 12 + TOpg 

, Jps*7*7* (g 2 , qf, o)Jps*7*7* (gi 1 gf 1 gig) „ nTc „ a j 

T „2 T , T1 2 t l i<TTp f ll t v\atpqiq2 

, ^PS'7'7* (gfi gl » ffl2)- 7r PS'7*7* (gfi gfi °) c „a„/3 r rc , 4 x 

H 2~T — 2 epAa/jgi q 2 e„apTq 2 , (34) 

<?2 m PS 

I 

where qi 2 = qi + q 2 - The function n( p ) p „ A(T is now only external photon momentum has been carried out. 
dependent upon two momenta since the limit k — >• for the 
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IV. RESULTS 

With the approach to hadronic LBL scattering within the 
functional approach outlined, and our truncation scheme de- 
fined we proceed to combine our propagators, vertices, ampli- 
tudes and form-factors together and calculate the respective 
contributions to hadronic LBL scattering in the muon g — 2. 



compatible with previous ones JH 0, which for the 
pion pole contribution are displayed in Table [TX| This is not 
surprising, since the form-factors themselves are compatible 
at a qualitative level, cf. Fig. [16] and all approaches make 
the assumption that pseudoscalar pole-dominance is valid far 
from the meson mass-shell. 



A. Pion-pole contribution to LBL 



B. Quark loop contribution to LBL 



To demonstrate parity between our approach and others in 
the determination of hadronic LBL, we calculate the ladder- 
exchange diagram of Fig. [3] assuming pseudoscalar pole- 
dominance. Once more, we re-iterate that on mass-shell this 
is identical to the pseudoscalar exchange diagram portrayed 
in Fig. |2l whereas off-shell we make the common assumption 
that the meson-exchange picture provides a good approxima- 
tion. 

In order to determine the pseudoscalar exchange contribu- 
tion we must numerically determine the dressed quark propa- 
gator, the quark-photon vertex, and the homogeneous Bethe- 
Salpeter amplitude for the pseudoscalar meson. Combining 
these together allows us to calculate from first principles the 
7T77 form-factor. We wish to emphasize again that the result- 
ing quark-photon vertex also contains time-like poles corre- 
sponding to vector meson exchange lf52tl . Thus the main ideas 
of VMD are naturally included here in the form-factor due 
to the non-perturbative approach that we employ. We have 
checked that the total numerical error of our calculation is of 
the order of one percent. In a similar fashion we also evalu- 
ate the corresponding form factors for the rj and ?/ mesons. 
We then use our results for the form factors to evaluate the 
pseudoscalar meson exchange contribution to LBL. For this, 
we use the off-shell prescription for the pseudoscalar Bethe- 
Salpeter amplitude proposed in Eq. d30l l for the exchanged 
mesons. This prescription gives a reduction of the contribu- 
tion that is similar to that found in other approaches. 

The systematic error of our calculation of the pseudoscalar 
exchange diagrams can be attributed entirely to the validity of 
the rainbow-ladder approximation, by the Maris-Tandy (MT) 
model, Eq. fTSY and the off-shell prescription for the mesons, 
Eq. d29l ). No other approximations have been used. While 
in the Goldstone-Boson sector the MT model works well, 
there is certainly a larger error in the flavour singlet sector. 
We therefore guesstimate a total systematic error: ten per- 
cent for the pion contribution, and twenty percent for the r\ 
and rj 1 contributions. With a numerical error of two percent 
we then obtain: ajjf L;7r ° = ( 57 5 ± g gj x 10 



-11 a LBL;i7 _ 



(13.6 ± 3.0) x 10~ n and a^ BL ^' = (9.6 ± 2.1) x 10" 11 
leading to 



,LBL;PS 



= (80.7 ± 12.0) x 10" 



(35) 



for the pseudoscalar meson exchange contribution to LBL. 
As compared to our previous work 12311 . the values for a^ BL;,) 

and a]f h ' v are slightly reduced due to a more consistent 
off-shell prescription in these channels. Our result (1351 ) is 



Having convinced ourselves that the method works we now 
focus on the quark loop contribution to LBL. Here we follow 
the strategy described around Eq. (O where the Ward identity 
obeyed by the four-point function with respect to the external 
field is exploited to construct quantities that are explicitly fi- 
nite. As for our numerical error we verified that we reproduce 
the well known perturbative result for the corresponding elec- 
tron loop with an accuracy of better than one per mille. In the 
quark loop we use the fully dressed quark propagators for the 
up, down, strange and charm quarks, extracted from the DSE, 
Fig. [TO] As already mentioned above, due to numerical com- 
plexity we are unfortunately not yet in a position where we can 
use the full, numerically determined quark-photon vertex used 
for the meson exchange contributions in the previous section. 
Instead, we use three different approximations to the full ver- 
tex and compare the results. As explained in section |HIB I due 
to the Ward-identity we are in possession of exact expressions 
for the non-transversal, Ball-Chiu part of the vertex, Eq. ( |2TI) . 
We exploit this knowledge to compare results with (a) a bare 
vertex, (b) the first term of Eq. (fUJi ( IBC) 



r M (p,<?) 



A(p2)+A(q2) 



7m 



(36) 



where p, q are the quark and antiquark momenta, and (c) the 
full Ball-Chiu (BC) expression Eq. (fJB. A comparison be- 
tween these three approximations may serve as a guide for the 
systematic error due to the relevance of vertex effects. We em- 
phasize, however, that only our most elaborate approximation, 
(c), satisfies the constraints of gauge invariance. Previous ap- 
proximations based on purely transverse parts of the vertex 
IU6I1 do not satisfy this constraint. We believe that the ansatz 
(c) provides an excellent basis for the calculation of the quark- 
loop diagram, which can and should be expanded in future 
work to also include transverse parts of the vertex. 
As a result of our calculation we find 



a 



a 



a 



LBL;quarkloop (bare vertex) 
LBL;quarkloop(lBC) 
LBL;quarkloop (BC) 



( 61 ±2) x 10" 11 
(107 ±2) x 10" 11 
(176 ±4) x 10~ n 



(37) 



for the quark loop contribution. Clearly these are sizable con- 
tributions. Whereas the bare vertex result roughly agrees with 
the number 60 x 10~ n given in fTill . the dressing effects of 
the vertex lead to a drastic increase. As compared to our result 
for the first part of the Ball-Chiu vertex B23I1 we again find a 
drastic increase from 107 x 10~ n to 176 x 10 -11 when the 
other two terms of the full Ball-Chiu vertex are included. In 
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Group 


Model 


aJ; BL (^°-pole) 


a^ BL (quark-loop) 


Bijnens, Prades, Pallante f 161 


EN1L 


59(11) 


21(3) 


Hayakawa, Kinoshita [10], HK and Sanda [8, 9] 


HLS 


57(4) 


9.7(11.1) 


Knecht and Nyffeler |43J 


LMD+V 


58(10) 


- 


Melnikov and Vainshtein [14] 


LMD+V 


77(5) 


t 


Dorokhov and Broniowski f 1 81 


NL X QM 


65(2) 


t 


Nyffeler [5] 


LMD+V 


72(12) 


t 


Our Result 


DSE 


58(10) 


136(59) 



TABLE H: Results for the n -pole and quark-loop contribution (where appropriate) to hadronic light-by-light scattering, in different models. 
For f the quark-loop correction is incorporated as a boundary condition on the pion-pole contribution, whilst for t. the quark-loop corrections 
are currently under investigation I58fl . 



this calculation we included effects from four quark flavours 
in the quark-loop. Their individual contributions are given by 



a 



a 



a 



LBUquarkloop (BC);u/d 
I" 

LBL;quarkloop (BC);s 
A" 

LBL;quarkloop (BC);c 



(158 ±3) X 10 

( 6 ± 1) x 10- 
( 12 ± 1) x 10 



-11 



-11 



(38) 



It is interesting to note that due to charge effects the heavy 
charm quark contributes more than the much lighter strange 
quark. 

We have checked the model dependence of the above re- 
sult by comparing with a similar calculation using a different 
model for the quark-gluon interaction IkToTl . The results are 
similar to the one in Eq. ( 1371 ) within an error margin of five 
to ten percent. Details will be given elsewhere. Due to these 
results we estimate an additional systematic error for our BC- 
result of 15 x 10~ n , which has to be added to the 4 x 10~ n 
given in Eq. ( 137b . 

In general, these large dressing effects also make it very 
hard if not impossible to guess the effect of the total vertex 
dressing without an explicit calculation. Certainly, however, 
given these findings, all previous estimates for the systematic 
error in the quark loop contributions seem to be an order of 
magnitude too small. 



V. CONCLUSIONS 

In this paper we have presented a new approach towards 
the anomalous magnetic moment of the muon. We have used 
a combination of Dyson-Schwinger and Bethe-Salpeter equa- 
tions to evaluate the pseudoscalar meson exchange contribu- 
tion and the quark loop contribution to LBL. Our only input 
is the Maris-Tandy model, a phenomenologically successful 
ansatz for the combined strength of the gluon propagator and 
the quark-gluon vertex. Our treatment of the meson exchange 
contribution to LBL is different from earlier approaches in 
that we do not rely on an ansatz for the PS77 form factor, 
but calculate this quantity starting from the basic equations 
of motion of QCD. Nevertheless, our result basically agrees 
with those from previous approaches. This result once more 
emphasizes that the meson exchange contributions to LBL are 
largely controlled by analytic constraints from QCD at large 
and small Q 2 . 



As for the quark-loop contribution, analytic constraints 
have been used which arise from the requirement of gauge 
invariance: the quark-photon vertex appearing in this loop has 
to satisfy the vector Ward-Takahashy identity. In contrast to 
previous approaches, we have implemented this identity by 
using the Ball-Chiu ansatz for this vertex. We believe this is 
a systematic improvement. The consequences are drastic: we 
observe a dramatic increase for the quark-loop contribution to 
LBL. Our result of (176±4) x 10~ u is more than three times 
larger than the constituent quark result of Ref. 11411 . 

When combining our two results, Eq. ( f35l > and Eq. d37i i, we 
arrive at a hadronic LBL contribution of 



a 



LBL;PS+quarkloop 



(257 ±31) x 10 



-11 



(39) 



This value, however, does not yet account for transverse parts 
of the quark-photon vertex in the quark-loop contribution 
and for effects from the right hand diagrams of Fig. [2] or 
Fig. [3] In general, it is difficult to gauge the effects of ad- 
ditional transverse vertex contributions in the quark-loop di- 
agram. In Ref. [[lfl part of these effects have been taken 
into account by using an ansatz motivated by vector-meson 
dominance (VMD) ideas. They found a reduction due to 
these effects of roughly 40 x 10 -11 . Since we agree with 
Ref. lfl6ll on the size of the pion exchange contribution, where 
VMD works very well, it may be justified to use their re- 
sult as a rough estimate for the size and also for the poten- 
tial error in these effects. We therefore add a contribution 4 

(-40 ± 40) x lO" 11 to arrive at 
: (136 ±59) x 10" 11 . 



of a 



LBL;quarkloop, transverse 



(I 



LBL;quarkloop (BC-t-transverse) 



The additional contributions due to the right hand diagrams 
of Fig. [2] or Fig. [3] are also difficult to judge. It may help, 
though, to observe that these involve an additional quark-loop. 
Typically such contributions are negative and of the order of 
ten to twenty percent of the leading- N c contributions 1I28H29I1 . 
Since on the other hand one also expects positive contributions 
of a similar size from non-pseudoscalar exchange diagrams 
Jit] we choose to subsume all these contributions to another 



Note that cross-terms between transverse and non-transverse contributions 
from the four vertices in the quark-box diagram could provide additional 
suppression. This is expressed in our error estimate. 
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a LBL;other = (q ± 20) x 10" 11 , where the error is clearly sub- 
jective. This gives us the following total hadronic LBL con- 
tribution 

aJf L = (217 ±91) x 1CT 11 , (40) 

in our approach. Note that the increase of the central value 
as compared to our previous result in Ref. ll23"ll is due to a 
combination of taking the full BC-vertex instead of 1BC and 
in addition accounting for the transverse corrections in the 
quark-loop using the results of Ref. fl6ll . Taken at face value 
these numbers together with the other contributions quoted in 
Ut] clearly reduce the discrepancy between theory and experi- 
ment. Combining our light-by-light scattering results with the 
other SM contributions gives: 

a theor. = H6 591891.0(105.0) x 10~ u . (41) 

To put this result in perspective we wish to recall the caveats 
that to our mind are tied to it. First, there is the contribution 
of transverse parts of the quark-photon vertex to the quark- 
loop diagram. Although the results of Ref. Il6ll may serve 



as an estimate, we definitely need to explicitly calculate these 
contributions in our approach. Second, there is the question 
whether the pseudoscalar meson exchange diagram provides 
for a good approximation of the gluon exchange contribution 
discussed around Fig. [3] Also this assumption needs to be 
questioned by an explicit calculation. In this sense, our results 
certainly do not provide final answers but still have to be seen 
as a further step towards a fundamental determination of . 

Finally, we point out that the current approach will also be 
checked by a calculation of the hadronic vacuum polarization 
contribution to a^. Preliminary results in this direction are 
encouraging. 
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